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Abstract. This paper uses the supernova data to explore the property of dark energy 
by some model independent methods. We first Taylor expand the scale factor a{t) 
and the luminosity distance to the fifth order to find out that the deceleration 
parameter go < 0. This result just invokes the Robertson- Walker metric. So the 
conclusion that the universe is expanding with acceleration is more general. Then we 
discuss several different parametrizations used in the literature. We also proposed two 
modified parametrizations. We find that cjdeo is less than —1 almost at la level from 
all the parametrizations used in this paper. We also find that the transition redshift 
from deceleration phase to acceleration phase is zt ~ 0.3. 
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1. Introduction 

The type la supernova (SN la) observation suggests that dark energy contributes 2/3 to 
the critical density of the present universe PQI3E]- SN la observation also provides the 
evidence of a decelerated universe in the recent past with the transition redshift z^ ~ 0.5 
01 El Ej- The cosmic microwave background (CMB) observations favor a spatially flat 
universe as predicted by inflationary models [3E]- There are many dark energy models 
proposed in the literature. For a review of dark energy models, see, for example, [H] 
and [lOj and references therein. However, the nature of dark energy is still unknown. It 
is not practical to test every single dark energy model by using the observational data. 
Therefore, a model independent probe of dark energy is one of the best ways to study 
the nature of dark energy. 

The type la supernovae (SNe la) as standard candles are used to measure the 
luminosity distance-redshift relationship d-^(z). So we can model the luminosity distance 
g?l to study the property of dark energy. Melchiorri etal. first found that dark energy 
may be a phantom type by combining different observational data to probe the behaviour 
of dark energy pi] . Huterer and Turner modelled the luminosity distance by a simple 
power law di,(z) = J2iLi c i zi [12] • Saini etal. used a more complicated function to 
model the luminosity distance ^3]. Another way to probe the nature of dark energy 
is to parameterize the dark energy equation of state parameter c^de- The simplest 
parametrization is the constant equation of state model u DE = constant. Several authors 
modelled cube as lube = Y^f=o^i z% [HI H3 EBJ • Apparently this parametrization is not 
good for high z. Recently, a stable parametrization o;de = + ^o^/(l + z ) was use d 
in [13 HH1 EE HI]. By fitting the model to SN la data, we find that u + uo a > 0, 
so this parametrization is not good at high z too. Jassal, Bagla and Padmanabhan 
modified this parametrization as ube — + w a z/(\ + z) 2 and the problem was solved 
because ujbe = at present and at high z [21] ■ More complicated functional forms for 
(jj-de( z ) were also proposed in the literature [23 123 123 1213 • We can also model the dark 
energy density itself. For example, a simple power law expansion VLbe = J2iLo Mz 1 was 
used to investigate the nature of dark energy [23 123 123 123 E3 EI] • There are other 
parametrizations, like the piecewise constant parametrization [33 E3 EH EH] • 

This paper is organized as follows. In section II, We first use a Taylor expansion 
to expand the scale factor, then we fit the model to the whole 157 gold sample of SNe 
la compiled by Riess etal. in p. By expanding the scale factor, the fitting parameters 
have physical meanings. In section III, we analyze the dark energy parametrization 
proposed by Alam etal. j2E]. In section IV, we first study the parametrization 
o>de = + u a z/(l + z) and point out that this parametrization is not good at high 
z. Then we study the parametrization ujbe = ojq + u a z/(l + z) 2 . In section V, we 
first investigate the parametrization proposed by Wetterich [2S], then we propose two 
modified parametrizations. In section VI, we give some discussions. 
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2. Model Independent Method 



In a homogeneous and isotropic universe, the Friedmann-Robertson- Walker (FRW) 
space-time metric is 

dr 2 



-dt 2 + a 2 (t) 



For a null geodesic, we have 



*o dt 



i a(t) 



n 



1 — kr 2 
dr 



+ r z dQ 



= fin), 
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(2) 



lo y/l — kr 2 

where 

sin -1 ri, k = 1 

ri, k = 

sinh -1 r 1; k = — 1. 

From Eq. (J2J), we get the luminosity distance d^ = a (l + z)r x by Taylor expansion 



(3) 



find 



0«L 



-(1 - q )z 2 + -( 9o + 3g 2 - 1 - Jo - ^fc)^ + 777(2 - 2<? 
2 o 24 



15g 2 - 15g 3 + 5 Jo + Wqaja + s + 2Q k + Qn k q )z 4 + 0(z 5 ), (4) 

where the redshift z is defined as 1 + z = ao/a(t), the subscript means that a variable 
is evaluated at the present time, the Hubble parameter H(t), the deceleration parameter 
q(t), the jerk parameter j(t) and the snap parameter s(t) are defined as 



1 da 
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(8) 



and flk = ^/( a o-^o)- The use of jerk parameter is equivalent to the statefinder used in 
|3*7| EHI • We may use the above expression (J1J) to probe the geometry of the Universe 
|39| liT)] . Note that the Taylor expansion of g?l may break down at high z and the 
actual behaviour of c?l rnay not be represented by finite number of terms. It is also 
straightforward to get 

H 2 (z) = H 2 [l + 2(l + q )z + (l + 2q +j )z 2 -^( So + q Jo)z 3 + 0(z%(9) 



q{z) = q -(q +2q 2 -j )z+ 4g 3 + 4g 2 + q - 2j, 



so 7j q 



Z 2 +O{z 3 ).{10) 



2 2 

Now let us find out qo, jo and Sq from the SN la data compiled by Riess etal., These 
parameters are determined by minimizing 

[/iobsOi) - V(Zi)} 2 



X 



E 



(11) 
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where (7j is the total uncertainty in the SN la observation and the extinction-corrected 
distance modulus n(z) = 5 log 10 (c?L(2 ; )/Mpc) + 25. In the fitting process, we use the 
SN la gold sample only and we consider a flat universe with = 0. Because we use 
Taylor expansion to get the luminosity distance, this expansion may break down at 
high z. Therefore, we first use the full 157 gold sample SNe, then we use those 148 SNe 
with z < 1.0. The best fit parameters to the whole 157 gold sample SNe are (go, jo, 
sq)=(— 1.1, 6.4, 39.5) with \ 2 — 174.2. The best fit parameters to the 148 gold sample 
SNe with z < 1.0 are (g , j , s )=(-1.7, 14.4, 149.4) with X 2 = 160.8. 

If we expand the luminosity distance d^ to the third order only, i.e., we only consider 
the parameters go an d jo in Eq. P|). then we find that the best fit parameters to 
the whole 157 gold sample SNe are: g = -0.64^11, jo = 1.2±\5 and \ 2 = 176.1. 
At 99.5% confidence level, g = — 0. 64^59, so we conclude that the expansion of 
the Universe is accelerating with 99.5% confidence. From Eq. ([10jh we get z^ = 
Qo/(<lo + 2gp — jo) = 0.595lo!i77- The best fit parameters to the 148 gold sample 
SNe with z < 1.0 are: g = -1.0 ± 0.4, j = 4.7j£J and x 2 = 161-3. At 99.5% 
confidence level, go = — 1.0t°;o, so we conclude again that the expansion of the Universe 
is accelerating with 99.5% confidence. With the best fit parameters, we find that 
zt = go/(go + 2go — jo) = 0.295lo;c)56. The contour plot for g and jo is shown in 
Figs. [UandEl 

So far our analysis uses the FRW metric only, we have not specified any gravitational 
theory yet. The above results are applicable to a wide range of theories. For example, 
g = (fi m o — 2^a)/2 and jo = fi m o + for the A-CDM model. If we expand the 
luminosity distance to the fifth order with the crackle parameter c(t) = (aH 5 )~ 1 d 5 a/dt 5 , 
then we need to add to Eq. (j3J) the following correction 

- ^(6 + 14g - 61g 2 - 160g 3 - 105g 4 + 110g jo 

+ 105g 2 j + 15g s + 27j - 10j 2 + Hs + Cq)z 5 . (12) 
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Figure 2. The plot of go and jo contour fitting to the 148 gold sample SNe with 
z < 1.0. 

The best fit parameters to the whole 157 gold sample SNe are (g , jo, s o> c o)=( — 1-1, 
7.6, 55.6, 676.6) with x 2 — 173.4. The correction to H^d-^ at z = 1.5 is about 
— 1.3 which is around 34%. The best fit parameters to the 148 gold sample SNe 
with z < 1.0 are (g , j , s , c ) = (-1.5, 11.4, 101.2,1484.2) with X 2 = 160.8. The 
correction to H^d-^ at z = 1.5 is about 0.19 which is around 8.5%. Therefore the 
introduction of the fifth order correction changes the value of go a little. We still have 
go < 0. It is clear that the kinematic determination of the cosmo logical parameters 
is better suited for low redshift SNe la. However, from the observational data, 
(z, fi(z)) = (1.4, 45.09), (z, n(z)) = (1.551, 45.3) and (z, fi(z)) = (1.755, 45.53), 
we see that A/i(z) = 0.21 when Az = 0.151 and A//(z) = 0.23 when Az = 0.205. 
Theoretically, we know that dfi{z) = (5/ ln(10))(<i d^(z) / d\Xz)dz)dz. From Eq. (JH), we 
get dfi(z) = (5/ ln(10)) (n / z)dz if the luminosity distance is dominated by the higher 
term z n . Combining the above analysis, we find that n ~ 1. Therefore, in this case, the 
higher term may not be the dominant term. 

We conclude that g < with 99.5% confidence. In other words, we conclude that 
the Universe is expanding with acceleration. 

3. "Taylor expansion" of Dark Energy Density 

In this section, we parameterize the dark energy density as j2Hj 

Q BE (z) =A + A 1 (l + z) + A 2 (l + z) 2 , (13) 

where Q DE (z) = 8nGp DE (z)/ (3H$), Q m0 = 8nGp m0 /(3H$) and A Q = 1-Q m0 - A x - A 2 . 
This parametrization is equivalent to Eq. Q with Q mQ = — (s + gojo)/3- The 
relationship between wde and z is 

l + z A 1 + 2A 2 (1 + z) 

UDE ~ 3 A + A 1 (l + z) + A 2 (l + z) 2 • 
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Figure 3. The best fit to the 157 gold sample SNe la with the prior fi m0 = 0.3 ± 0.04. 
The upper panel shows cjde(z), the dotted dash lines are the la regions. The lower 
panel shows Q m (z) and Qde(z) 

With the above parameteriaztion, we find that f2 DE <C Q m and u DE « —1/3 when 
z > 1. Combining the above two equations, we find that the transition redshift z^ 
satisfies 

Q m0 (l + zf -A x {l + z)- 2A = 0. (14) 

The best fit to the whole 157 gold sample SNe gives A x = —5.79, A 2 = 2.9 and fl m0 ~ 
with x 2 — 174.05. If we use a Gaussian prior Q m0 = 0.3 ± 0.04 |4~T] . then we get the 
best fit parameters A\ = —4. 21^4 and A 2 = l-7liJ with x 2 = 174.21. Substitute these 
parameter values to Eq. (fT4|) . we find that zt = 0.35. The evolutions of the dark energy 
density and c^de are shown in Fig. El Alam et al. showed that the SNe la data favored 
an evolving dark energy model by using the above reconstruction HZj • They also 
showed that zt ~ 0.4. Our results are consistent with those analysis. 

Because it is possible that c^de < — 1 ; so we consider another two parameter 
representation of dark energy 

n DE (z) =B + B^l + z)+ B_i/(1 + z), (15) 

where B = 1 — Q m0 — B\ — B_\. with this parametrization, we get 

1 B^l + zf-B^ 

^ DE 3B l (l + z)* + B (l + z)+B_ 1 ' 

The above equation tells us that Ode <S and u DE rs — 2/3 when z ^> 1. Combining 
the above two equations, we find that the transition redshift z^ satisfies 

Q wa (l + z) 3 -B 1 (l + z)-2B -^^ = 0. (16) 

1 ~r Z 

The best fit to the whole 157 gold sample SNe la gives B_i = 6.87, B\ = 6.14 and 
^mo ~ with x 2 — 173.2. If we use a Gaussian prior Q m0 = 0.3 ± 0.04, then we get the 
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Figure 4. The best fit to the 157 gold sample SNe la with the prior £) m0 = 0.3 ± 0.04. 
The upper panel shows cjde(z), the dotted dash lines are the la regions. The lower 
panel shows Q m (z) and £Id&(z) 



best fit parameters £>_i = ^ll^ and B\ = 2. 8^0 with x 2 — 173.65. Substitute the 
best fit parameters into Eq. (fTtjj) . we get z T = 0.30. The evolutions of u DE and Q-^e are 
shown in Fig. 0] 



4. Stable Parametrization 



In this section, we first consider the parametrization [TT| ITS] 



^DE — 1^0 + 



l + z' 



;i7) 



When z ^> 1, we have c^de ~ + <-<V The dark energy density is 

^de = ^deo(1 + zf( 1+ " 0+ ^ exp(-3u a z/(l + z)). 



Combining the above two equations, we find that zt satisfies 



H m o + ( 1 — 



mOj 



1 + 3CJ + 



l + Z 



:i+z 



exp 



SuJ n Z 



0.(18) 



l + z 

The best fit to the whole 157 gold sample SNe la gives ojq = —2.5, u a = 3.7 and 
^mo = 0.46 with x 2 — 173.5. If we use a Gaussian prior Q mQ = 0.3 ± 0.04, then we get 
the best fit parameters ujq = — 1.6loJ an d u a = 3.3l|; 7 with x 2 — 173.92. Substitute 
the best fit parameters into Eq. (fTSjl. we get zt = 0.35. The evolutions of ube and Qbe 
are shown in Fig. 

From Fig. we see that the dark energy density is greater than the matter density 
at high z because u + uo a > 0. So this stable parametrization may not be a good 
choice at high z. Recently, Jassal, Bagla and Padmanabhan considered the following 
parametrization [2~T] . 

LU„.Z 



ujde — ^0 + 



'l + z 



,2 ' 



(19) 
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Figure 5. The best fit to the 157 gold sample SNe la with the prior ft m0 = 0.3 ± 0.04. 
The upper panel shows u>x>e(z), the dotted dash lines are the la regions. The lower 
panel shows Q m {z) and Qde(z) 



When z ^> 1, we have wde ~ Hi- The dark energy density is 

^de = ^deo(1 + z) 3(1+ " o) exp(3^ 2 /2(l + zf). 
Combining the above two equations, we find that zt satisfies 

n m0 + (i-n m0 ) ^i + 3c + ^^ (i + ^exp { ^f z y j =o-(2Q) 

The best fit to the whole 157 gold sample SNe la gives ojq = —2.5, u a = 7.6 and 
Q m o = 0.42 with x 2 — 173.3. If we use a Gaussian prior fl m0 = 0.3 ± 0.04, then we get 
the best fit parameters cu = — 1.9t?j and cu a = 6.6 ± 6.7 with \ 2 — 173.41. Substitute 
the best fit parameters into Eq. (|20|). we get z^ = 0.30. The evolutions of u DE and 
Qbe are shown in Fig. El From Fig. |f3 it is clear that the dark energy density did not 
dominate over the matter energy density at high z. Our result is consistent with that 
obtained in 



5. Wetterich's Parametrization 

In this section, we first consider the parametrization given in 

ude = [i+b\T(i+z)r (21) 

When z 3> 1, we have urjE ~ 0. The dark energy density is 

n DE = n DE0 (i + z) 3+3 ^ 1+bl ^l 

Combining the above two equations, we find that z^ satisfies 

+ (1 - IU) (l + |1 + M '( 1 ° + 2)|2 ) (1 + 2 )W| 1 +m»< 1 +,>1 = o. (22) 
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Figure 6. The best fit to the 157 gold sample SNe la with the prior fi m0 = 0.3 ± 0.04. 
The upper panel shows loue(z), the dotted dash lines are the la regions. The lower 
panel shows Q m (z) and flw(z) 



-1 _ . 

2 /I 



0.2 0.4 0.6 0.8 



1 1.2 1.4 1.6 1.8 




Redshift z 



Figure 7. The best fit to the 157 gold sample SNe la with the prior J7 m0 = 0.3 ± 0.04. 
The upper panel shows loy>^{z), the dotted dash lines are the la regions. The lower 
panel shows Q m (z) and £Ide(z) 



The best fit to the whole 157 gold sample SNe la gives uj = —1.84, b = 5.85 and 
^mo ~ with x 2 — 173.09. If we use a Gaussian prior Q m0 = 0.3 ± 0.04, then we get 
the best fit parameters cu = — 2.5t\'^ and b = 4.0l3 1 5 4 with x 2 = 173.15. Substitute the 
best fit parameters into Eq. (j22J), we get = 0.26. The evolutions of ojde and ^de are 
shown in Fig. [7| Because the best fit of the above parametrization gives Q m o ~ which 
is not physical, we first modify the above parametrization as 

UJ 



^DE 



1 + 61n(l + z) 



(23) 



Model independent analysis of dark energy: Supernova fitting result 



10 



co DE (z)=co /(1+bln(1+z)) 



3 -3 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 




Redshift z 



Figure 8. The best fit to the 157 gold sample SNe la with the prior fi m0 = 0.3 ± 0.04. 
The upper panel shows ojde(-z), the dotted dash lines are the la regions. The lower 
panel shows fl m (z) and Qde(z) 



When 2 > 1, we have u DE ~ 0. The dark energy density is 

fi DE = ^deo(1 + ^) 3 [1 + 61n(l + z)f^ /b . 
Combining the above two equations, we find that z^ satisfies 



^mo + (1 — Q 



mOJ 



1 



[1 + b\n(l + z)] c 



0. 



(24) 



l + Mn(l + z) 

The best fit to the whole 157 gold sample SNe la gives u>o = —3.05, b = 36.8 and 
^mo ~ with x 2 — 172.75. If we use a Gaussian prior Q m0 = 0.3 ± 0.04, then we get 
the best fit parameters u = -3.4i?£ 7 and u a = 17.8lijj?4 3 with x 2 = 172.91. This 
modification does not solve the problem of f2 m o ~ 0. Substitute the best fit parameters 
into Eq. pijl . we get z^ = 0.25. The evolutions of cjde and f^DE are shown in Fig. El 
Now let us consider another modification 



^DE — + 



l + ln(l + z)' 

When z ^> 1, we have wde ~ ^o- The dark energy density is 



(25) 



n 



DE 



n 



DEO 



l + Z 



i3(l+wo) 



[1 + ln(l + z)] 



3u> a 



Combining the above two equations, we find that zt satisfies 

3u; a 



^m0 + (l — ^ 



mOJ 



1 + 3^ 







l+^) 3 ^°[l+ln(l+2 



l3a; a 



0.(26) 



1 + Mn(l + z) 

The best fit to the whole 157 gold sample SNe la gives u = 2.2, oj a = —4.7 and 
O m o = 0.454 with x 2 — 173.47. If we use a Gaussian prior Q m0 = 0.3 ± 0.04, then we 
get the best fit parameters cu = 2.4t|g and cu a = — 4.111;! with x 2 = 173.81. Substitute 
the best fit parameters into Eq. (J2EJ), we get z^ = 0.34. The evolutions of o>de and ^de 
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Figure 9. The best fit to the 157 gold sample SNe la with the prior Q m0 = 0.3 ± 0.04. 
The upper panel shows lo^{z), the dotted dash lines are the la regions. The lower 
panel shows fl m (z) and fljm(z) 

are shown in Fig. Although this modification solves the problem of Q m o ~ 0, it is not 
good at early times because the dark energy density dominated over the matter energy 
density at early times as shown in Fig. El 

6. Discussions 

The SN la data shows that the expansion of the Universe is accelerating. This conclusion 
derived from Eqs. (HJ) and (fT2|) does not dependent on any particular model. We used 
the parametrizations (fT3|) . (fTB]) . (|T7|) . (fT^)l and (^1]) proposed in the literature to discuss 
the property of dark energy. We also proposed two modified parametrizations (J23|) and 
(OB")) . By using the above parametrizations, we derived the equations satisfied by the 
transition redshift. In order to see the property of ude(z), we re-plot ude(z) for all the 
models considered in this paper together in Fig. E3 From Fig. we see that: (a) 
<^deo < — 1- This is also true at la level. So the current SN la data seems to marginally 
favor the dark energy metamorphosis suggested in [2E1 EZl • This does not mean that we 
can exclude the A-CDM model; (b) uy)e(z) increases when z increases. ujt>v{z) changes 
more rapidly at low z than at high z. This property may be due to the choice of the 
parametrizations we made; (c) ~ 0.3. We also see that the parametrization (fTT?|) is a 
good choice. It avoids the problem that the dark energy dominated the matter energy 
at early times and the best fit Q m o to the SN la data for this parametrization is not 
close to zero. The problem of f2 m0 ~ is not a serious problem because x 2 depends 
on Q m0 weakly for all the models discussed in this paper. Daly and Djorgovski found 
that z T ~ 0.4 by using a model independent analysis [231211 • In our analysis, we used 
Friedmann equation and some priors to interpret the SN la data. As shown in j!2] , the 
interpretation of the observational data changes drastically if the priors are removed. 
We would like to stress that the results obtained in this paper are consistent with other 
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Figure 10. The evolution of ude for different parametrizations. 
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